Astrophysical searches for a hidden-photon signal in the radio regime 
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Common extensions of the Standard Model of particle physics predict the existence of a "hidden" 
sector that comprises particles with a vanishing or very weak coupling to particles of the Standard 
Model (visible sector). For very light (m < 1CP 14 eV) hidden U(l) gauge bosons (hidden photons), 
broad-band radio spectra of compact radio sources could be modified due to weak kinetic mixing 
with radio photons. Here, search methods are developed and their sensitivity discussed, with specific 
emphasis on the effect of the coherence length of the signal, instrumental bandwidth, and spectral 
resolution. We conclude that radio observations in the frequency range of 0.03-1400 GHz probe 
kinetic mixing of ~ 10 -3 of hidden photons with masses down to ~ 10 _ 17 eV. Prospects for improving 
the sensitivity with future radio astronomical facilities as well as by stacking data from multiple 
objects are discussed. 



I. INTRODUCTION 

Finding experimental evidence for physics beyond the 
Standard Model (SM) of particle physics is one of the pin- 
nacles of present-day physical research, embracing both 
extensive laboratory studies and indirect (primarily) as- 
trophysical measurements made across a very broad 
range of energies. Most of present-day SM extensions 
into a more generic, unified scenario predict existence of a 
class of particles only weakly interacting with the normal 
matter: weakly interacting massive particles (WIMP) at 
a mass scale of 0(100) GeV [J41 and ultralight weakly 
interacting sub-eV particles (WISP) [HI]. 

The existence of ultralight particles has been argued to 
be at least theoretically plausible in a number of different 
scenarios including additional pseudo-scalar (axions and 
axion-like particles (ALP), </>, |9jjjdJ_) as well as vector 
fields (U(l) hidden photons, 7 S , [Hi, H ) • In either 
scenario, the prevailing non-baryonic matter could be ex- 
plained by these ultra-light fields [H, [l4| . 

Hidden photons arise in low-energy extensions of the 
SM which leave all the SM fields uncharged under the ad- 
ditional hidden gauge group. Interaction between hidden 
photons and massive SM particles is expected to be sup- 
pressed by the particle masses. However, it can be mani- 
fested by kinetic mixing with nor mal photons 0, [II III. 
In low-energy SM extensions with hidden photons, the 
kinetic mixing is expressed by the effective Lagrangian 
describing two-photon interactions L = £sm + Ai + 
[cf., [HI [HI [lTj, with £sm denoting the Maxwell La- 
grangian for the SM photon field, the £h term describing 
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the Proca-Lagrangian for the hidden photon field, and 
C x representing a gauge-invariant kinetic mixing term. 
The kinetic mixing term induces photon oscillations be- 
tween the massless "normal" state (7) and a non-zero 
mass "hidden" state (7s). In this hidden state, photons 
acquire a non- vanishing mass and propagate on space-like 
geodesies, without any interaction with normal matter. 
The physical properties of a hidden photon can be com- 
pletely described by its mass m 7 „ and the kinetic mix- 
ing with an SM photon (expressed by the mixing angle 
x) ■ Theoretical predictions for y are falling in the broad 
range between 10~ 16 and 10~ 2 [3, LH423|- 

Accelerator experiments are generally optimized to 
search for new heavy particles such as WIMPs, and 
therefore of limited sensitivity and mass reach for 
WISPs. Hence, the potential discovery of ultralight par- 
ticles requires high-precision experiments for which non- 
accelerator setups often appear more promising [23H25j . 

Evidence for the 7~7 S oscillations has been searched 
for in a number of laboratory 26H3l( and astrophysi- 
cal experiments 0, El HI [HUfocusing inp articular 
on "light shining thr oug h the wall" (LSW, J36|) experi- 
ments such as ALPS [27| and searches for hidden photons 
from the Sun (SHIPS; [H, H3). The non-detection of 
hidden-photon signals has so far yielded strong bounds 
on the kinetic mixing parameter \ f° r a broad range of 
hidden photon masses @, 0, E3, and references therein] . 
The mass range which has been currently probed extends 
down to m 7a = 2 x 10 _14 eV, with the lowest hidden- 
photon masses probed by the WMAP CMB measure- 
ments in the radio domain at frequencies above 22 GHz 
[H[. Below 10" 14 eV, only weak limits of x ~ 10" 2 H 
have been obtained from analysis of early measurements 
of magnetic field around Earth and Jupiter [H| , and no 
limits are reported for m 7a < 5 x 10 -16 eV. 

Radio observations at frequencies below 22 GHz offer 
an excellent (if not unique) tool for placing bounds on 
the mixing angle x f° r m 7„ < 10 _14 cV. Initial bounds 
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on x can be obtained from existing radio data on com- 
pact, weakly variable objects with well-known radio spec- 
tra (such as young supernova remnants (SNR), plan- 
etary nebula, and steep spectrum radio sources typi- 
cally used for the absolute flux density calibration of 
radio telescopes). With this approach, one can reason- 
ably expect to reach \ < 0.01. Propagation through 
a refractive medium (in which an SM photon also ac- 
quires an effective mass m 7 ) can strongly affect this limit, 
improving it substantially near the resonance condition 
?7i 7a = m 7 , and sup pres sing the hidden photon conver- 
sion at m 7 3> 77i 7a [33J . The latter effect can become 
important at to 7s < 10 -15 eV. 

Placing better bounds on \ can now be achieved by 
using the expanded capabilities of existing radio tele- 
scopes (utilizing the upgraded broad-band coverage and 
spectral resolution of the Effelsberg 100-meter antenna 
and the Karl Jansky Very Large Array (JVLA) in the 
0.3-40 GHz frequency range), by extending the measure- 
ments both to lower frequencies (0.03-0.3 GHz) covered 
by the Low Frequency Array (LOFAR) and to submil- 
limeter wavelengths probed by Atacama Large Millimeter 
Array (ALMA), and by employing the superb brightness 
sensitivity (~ 1 ujy [39[) of the SKA[i(| and its precur- 
sors, MeerKAT[4l[ and ASKAP[42[. 

The existing data for the primary absolute flux density 
calibrators in the radio regime (such as Cas A, Tau A, and 
Cyg A) [43| feature a few dozens of absolute flux density 
measurements made in the 0.01-30 GHz spectral range 
and reaching a ~ 3-5% accuracy. These data should en- 
able placing a bound of x ~ 0.02, from measurements 
of the r.m.s. of deviations from a canonical source spec- 
trum. If the spectral resolution is sufficiently high to 
assess the periodicity in the oscillation signal (particu- 
larly at the lower end of the spectrum), both the limits 
on x an d the range of photon mass studied can be im- 
proved. Further improvements of the bounds on x can be 
achieved by stacking the signal from a number of objects, 
under the condition that the observations of different ob- 
jects are sensitive to the same range of the hidden-photon 
mass. 

In this paper, a methodology and prospects for detec- 
tion of the hidden photon signal in the radio regime are 
considered. The basic physics of the 7 — 7 S oscillation 
and the propagation of the hidden-photon signal are de- 
scribed in Section [II] Methods for the detection of the 
oscillation signal with different instruments and targets 
are discussed in Section Hill and potentials of these studies 
are discussed in Section [TVl 



where A^ u and -B M „ are the respective field-strength ten- 
sors. The non-zero kinetic mixing angle x implies a mis- 
match between the interaction and propagation eigen- 
states, which induces oscillation between the two states 
(with a close analogy to the neutrino oscillation effect 
[44j ) . The last term accounts for massive hidden photons 
via Higgs or Stiickelberg mechanisms, where the former 
case suffers from additional constraints @ . The probabil- 
ity of the 7 — > 7 S conversion after propagating a distance 
L in vacuum is then given by 
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7->7=( L ) =a x sin \j^L 
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where all quantities are expressed in natural units, and E 
and v are the energy and frequency of the normal photon 
El El HI- It should be noted that the validity of Eq. [2] 
is restricted to the case of m 7s <C 2-kv, which is fulfilled 
for all considerations in this paper. 

The first term of Eq. [5] describes the amplitude of the 
oscillation, a x = sin 2 (2x) ~ 4% 2 (for ^Cl), that can 
be identified as a periodic signal over a range of distances 
L or wavelengths A = The second term, 
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gives the periodic signature of the oscillation. The os- 
cillation signal affects the broad-band spectrum, F(v), 
of an astrophysical source, which results in a received 
spectrum F 1b (v) = F(u) (1 — P 7 _> 7s ). 

The resulting spectrum F^ s (y) will have local minima 
and maxima at the frequencies v m i n ,i = f*/(2 i — 1) and 
"max,! = i/*/(2i), i e N, where 
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is the frequency of the first (highest frequency) minimum 
obtained with i = 1. The frequency, i/ ma x,i, of the first 
maximum defines the characteristic wavelength, 
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of the periodic modulation induced by the hidden-photon 
signal on a broad-band spectrum F(X) in the wavelength 
domain. 



A. Oscillation and coherence lengths 



II. PHOTON OSCILLATIONS IN THE RADIO 
REGIME 

For a photon field, A^, and a hidden photon field, B^, 
the kinetic mixing term is given by (lH | 



The oscillation length of the hidden-photon signal is 
set by the condition y> 7g (L, v) = tt, 
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Dccohercncc effects in the photon and hidden-photon 
mass eigenstates may arise during propagation, owing to 
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the finite mass, m 7a , of the latter. An accurate quan- m 7 , the oscillations are rapidly damped (medium sup- 



tum mechanical treatment of the oscillation probability 
P 7 _> 7b using wave packets (Eq. [2]has been derived assum- 
ing plane waves) yields an upper bound on the accessi- 
ble distance range [II El El, L coh = 4V2a x E 2 jm\, 

where a 2 = 2 + a® 2 denotes the quantum mechanical 
uncertainties of the production and detection processes, 
respectively. The non-thermal radio emission of compact 
radio sources is produced by synchrotron radiation, with 
the mean energy loss path of synchrotron emitting elec- 
trons = At (note that the mean energy loss path is 
of the same order of magnitude as the gyro radius of the 
relativistic electrons under consideration). A reasonable 
estimate of the cooling time At = 2kv c /(— dE/dt) can 
be obtained at the critical frequency v c of synchrotron ra- 
diation (averaging over the pitch angle) [47] . eventually 
yielding 
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where B denotes the magnetic field inside the considered 
source. To avoid a freeze out of vacuum 7 — 7 S oscilla- 
tions, the probed distance L therefore needs to fulfill the 
condition L osc < L < L C oh- 

Together with the condition i osc < L < L co h, Eqs. [6] 
and [7] imply effective lower and upper bounds on the 
hidden-photon mass that can be probed with radio data 
of a compact synchrotron emitting source at distance 
Lq above a given frequency vq ■ Table [I] lists corre- 
sponding bounds considering three distinct types of ra- 
dio sources, namely nearby Galactic supernova remnants 
(SNR), distant radio-emitting lobes of active galactic nu- 
clei (AGN), and AGN cores at cosmological distance 
scales. It demonstrates that the oscillation length and 
decoherence effects should enable effective radio mea- 
surements for a range of hidden-photon masses between 
-3 x l(T 19 eV and -4 x lfr n eV. 



B. Propagation through refractive media 



pression), with Xr — > tt/2 for £ 3> 1. The condition 
sin 2xr > sin 2\ implies that hidden photons with m 7a 
and x can be detected in a medium with the effective 
photon mass m 2 < 2 m 2 a cos 2%, which approximately 
yields m 7 < s/2m le for x <C 1. 

For photon propagation in the interstellar (ISM) and 
intergalactic (IGM) medium, the dominant factor is scat- 
tering off free electrons and neutral atoms (with the 
medium described by the electron and proton number 
densities, n a and n p ). In this case, the effective pho- 
ton mass, m 2 w uip — 2uj 2 {n — l) me dium [HI, depends 
on the photon frequency u, the plasma frequency Wp = 
4:iran e /m e , and the refraction index (n — l) m edium of the 
medium. 

For Galactic objects and extragalactic objects at small 
rcdshifts (z < 1), the medium can be assumed strongly 
ionized (with the ionized fraction of hydrogen X c = 
n e /n p — > 1). Contributions from helium and heavier ele- 
ments can be neglected. This yields [HI 
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For observations in the radio band, m 7 ~ wp ~ 3.7 x 
10~ n (nc/cm -3 ) 1 / 2 eV provides a good estimate of the 
effective photon mass under the assumption of X e sa 
1. This limits, formally, the hidden-photon mass that 
can be probed in Galactic and extragalactic objects to 
~ 10~ 13 eV and ~ 10~ 14 eV, respectively (for generic as- 
sumptions for the average densities of tiigm ~ 10 -5 cm" 
HI and niGM ~ 10~ 7 cm~ 3 (H). However, both the 
ISM and the IGM arc inhomogencous, with variations of 
the density exceeding 2-3 orders of magnitude [50|, H3 ■ 
Hence, generation and propagation of the hidden-photon 
signal depend on the line-of-sight (LOS) properties of the 
medium and its inhomogeneities. 



Photon propagation through a medium with refractive 
index n can be described by introducing an effective pho- 
ton mass m 7 to the Lagrangian C. This operation affects 
the kinetic mixing term, and the resulting effective mix- 
ing angle Yr depends on the mass ratio £ = m 2 /m 2 so 
that Ql M, HI 



sin 2xr 



sin2x 



sin 2 2x + (cos2 X -£) 2 



(8) 



For small effective photon masses, { « 1, the 7-7S os- 
cillations approach the vacuum regime, with Xr — > X- 
A resonance with the maximum amplitude Xr = 7r /4 
of the oscillations is reached at the resonant mass ra- 
tio £ = cos2x- For higher effective photon masses, 



C. Effects of inhomogeneous media 



In an inhomogeneous medium, the minimum value of 
m 7s without medium damping is limited by the elec- 
tron density of underdense regions with n c <C n c . The 
traversed underdense region needs to be sufficiently ex- 
tended to affect the propagation, such that the LOS path- 
length, Zi os , in these regions is ^ os 3> L osc (for observa- 
tions in the radio domain at v > 100 MHz, this effectively 
limits m 7s > 10 -16 eV for Galactic objects, while permit- 
ting probing hidden photons with m 7s > 10 _19 eV with 
extragalactic targets, see Table HJ. The resulting spectral 
pattern after propagating through underdense regions re- 
mains "frozen" during subsequent propagation through 
denser regions (where m 7s > 2.6 x 10 _11 (ri c /cm -3 ) 1 / 2 ), 
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TABLE I. Lower and upper bounds on the mass range of hidden photons, m 7a and m 7B , respectively, detectable with radio 
data (y > 1 GHz) of SNR and AGN. For m 7s < m' i; no 7~7 S oscillations will arise; for m 7a > m 7a , decoherence effects yield a 
freeze out of the oscillation signal. The columns B and AL list the assumed magnetic field and distance range, respectively. 





B 


AL 




m 7 






m 7a 


Source class 


[mG] 






[eV] 






[eV] 


SNR 


0.1 


1-10 kpc 


2 x 10" 


10 -6 x 10" 


-its 


1 x 10" 


11 -4 x 10" 11 


AGN: Radio lobes 


1 


0.02-3 Gpc 


3 x 10" 


19 -4 x 10" 


-18 


3 x 10" 


15 -3 x 10" 14 


AGN: Nuclear regions 


10 3 


0.02-3 Gpc 


3 x 10" 


19 -4 x 10" 


-18 


10" 


-16_ 10 -15 



as both direct and reverse photon conversions are sup- 
pressed there. Therefore, for Galactic objects, the low- 
est detectable hidden-photon mass would be achieved for 
photon beams propagating between the Galactic arms, 
while propagation through cosmic voids would set the 
lower limit on the hidden-photon mass that can be de- 
tected in the broad-band spectra of extragalactic targets. 

A density of free electrons n fi i nn ~ 0.005-0.01 cm -3 
is measured in the local ISM [52|, [53| , and there is am- 
ple evidence for n c to vary strongly across the Galaxy 
[52l [53| . with n c > 10 cm -3 near the Galactic center, 
n c r~j 10"" 2 cm"" 3 in the spiral arms, and n e <C 10 — cm 
above the Galactic disk. In mini "void" regions of ~ 1 kpc 
in extent and located between the spiral arms [53| , n c < 
10~ 6 cm" 3 can be found [5l|, which is similar to the val- 
ues typically measured in the IGM. Hence, detectability 
of hidden photon oscillation in Galactic sources should 
depend strongly on the LOS to a specific target, with 
likely m 7a > 10~ 12 eV detectable for LOS not crossing 
the spiral arms, while to 7s > 10~ 15 eV may still be de- 
tectable for objects at high galactic latitudes, and the 
LOS crossing inter-arm plasma and the local "voids" . 

The electron and proton densities in different struc- 
tural components of the IGM can be estimated from 
observations, i.e., IRAS data [54| and SDSS data [55| . 
as well as detailed numerical simulations of large-scale 
structures [H [H H3- The SDSS data indicates that 
the voids have a volume filling factor of 0.62 and a me- 
dian size of 17 h~ l Mpc, and Ref. [Ho[ find the bary- 
onic matter density, fib, void ~ (0.045 ± 0.015) fib, where 
h = iT/(100 kms" 1 Mpc -1 ) is the dimensionless Hubble 
constant and fib = 0.046 ± 0.002 is the average baryon 
density in the Universe [58| . Table |H] summarizes the 
mass fraction, m/rrib, the volume fraction, V/V c , and the 
relative density, p/pb, of different IGM components (the 
warm, warm-hot ionized (WHIM), hot, and void compo- 
nents, measured with respect to the total baryon mass 
TOb, density pb, and the comoving volume V c ). Based on 
these values, one can estimate the electron density in the 
individual IGM components 

_ y- Pc^b Pmcdium ( A C\ \ 

^c. medium — ^-c \ TT 3 V J 

TTip + X c m c p h 

where p c = 3H/(8irG) is the critical density of the Uni- 
verse, m p and m c are the proton and electron masses, 
and G denotes the gravitational constant. This yields 
n e (pb) = 2.5 x 10~ 7 cm"~ 3 and enables calculating n e 



and respective limits on m lB for different baryonic mat- 
ter components as listed in Table [TT] (compiled from the 
results reported in [5fl I54T457T ] ) . 

TABLE II. Baryonic matter components and hidden photon 
propagation. Each IGM component is described by the tem- 
perature T, mass fraction m/mb, comoving volume fraction 
V/V c , local density relative to the average baryon density 
p/pb, estimated average electron density n c , and resulting 
minimum detectable hidden-photon mass m 7s as estimated 
from the average electron density. 



Baryonic 


T 


m/mb 


V/Vc 


p/Pb 


n c 




m 7s 




component 


[K] 








[cm -3 " 




[eV] 




Galaxies 


<W 


0.054 


0.002 


27.0 


6.7x10" 


-6 


6.7x10" 


-14 


Warm IGM 


<10 5 


0.350 


0.342 


1.02 


2.6x10" 


-7 


1.3x10" 


-14 


WHIM IGM 


<10 6 


0.471 


0.030 


15.7 


3.9x10" 


-6 


5.1x10" 


-14 


Hot IGM 


>10 6 


0.097 


0.006 


16.1 


4.0x10" 


-6 


5.2x10" 


-14 


Voids f 


-10 6 (?) 


0.028 


0.620 


0.04 


1.1x10" 


-8 


2.7x10" 


-16 



In cosmic voids, strong density gradients are observed 
[55I l59l and references therein] , with densities at the void 
center being at least 2-3 orders of magnitude lower than 
at the edge of the void. Given the galaxy density as a 
tracer of the gas density, the density profile of electrons, 
n c (r' ) , can be calculated from the average radial density 
profile obtained from the galaxy counts (see Fig. 4-6 in 
|55|), where r' = r/r vo ;d, with r vo id denoting the void 
radius. The profile was normalized to reproduce the void 
average density reported in |50| . For such a profile, the 
minimum m 7a can then be calculated by evaluating n e (r) 
at Losc/2 (to account propagation on both sides from the 
center of the void). The resulting values are given in 
Table IIIII for several typical values of the void size and 
observing frequencies. 

This calculation demonstrates that for the assump- 
tions used here, hidden photons with masses as low as 
~ 10~ eV can in principle be probed with astrophys- 
ical measurements of targets located behind sufficiently 
large voids. A more detailed account of ISM/IGM inho- 
mogeneities on generation and propagation of the hidden 
photon signal would rely on extensive numerical simula- 
tions of large scale structures, which is beyond the scope 
of this paper. In the following discussion, the generic 
to 7s limits obtained for Galactic (> 10~ 15 eV) and ex- 
tragalactic (> 10 -17 eV) objects will be assumed, while 
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TABLE III. Minimum detectable hidden-photon mass m 7a in 
eV, assuming photon propagation through cosmic voids. 



Void 


Observing frequency 




2r void [Mpc] 


30 MHz 100 MHz 1 GHz 




10 


1.7x10"" 2.3xlO" iY 4.1x10" 


17 


30 


1.3xl0" 17 1.7xl0" 17 3.1x10" 


17 


100 


9.6xl0" 18 1.3xl0" 17 2.3x10" 


17 



In the following discussion, it is assumed that the ob- 
servational setup for a broad-band spectrum measure- 
ment can be simplified and characterized by a range 
\y\ ; ^2] of frequencies probed with a spectral resolution 
Av (implying that the flux density measurements are 
made at average intervals Av). In general, Av may vary 
across the frequency range, hence it is used here only in 
the sense of defining the total number of independent flux 
density measurements N mcs = (V2 — v\)j Av. 



performing calculations for the m 7 < y/2m^ e case. 

It is interesting to note that the signal from hidden 
photons with a sufficiently low mass, m-^ ee , should be 
unaffected by propagation through a high-density struc- 
ture with a characteristic size Z s , for which l s < L osc . This 
implies m^ 66 = 2ny/2 (v/l s ) 1/2 and yields, at v = 1 GHz, 
m^ s ee « 1.3 x 10" 16 
(l s w 20kpc). 



eV for propagation through galaxies 



III. DETECTION OF THE OSCILLATION 
SIGNAL 

For measurements in the m 7 < \[2ra ls regime, the two 
main factors limiting the sensitivity for a hidden photon 
signal are the spectral energy distribution of the (typi- 
cally, multicomponent) astrophysical signal and spectral 
range covered by the resolution of astronomical instru- 
ments. As the hidden-photon signal modulates the as- 
trophysical signal, the latter has to be well understood 
before attempting to detect the oscillation signal in a 
broad-band spectrum. The astrophysical signal can be 
modeled with JA(v), such that a condition 



KS V ) = F n«(y)/M{v) = c M (i- p n ) + 



(ii) 



is achieved (or approached), where Cm is a constant (ex- 
pecting Cm ~ * 1) and cr rm s is the residual fractional noise 
due to measurement errors and systematic uncertainties 
of the fit by M(v) (with er rms <C Cm)- In the radio 
regime, measurement errors will be dominated by the 
system noise and atmospheric/ionospheric fluctuations, 
while the effect of scattering in interstellar and inter- 
galactic plasma would be orders of magnitude smaller 
and could be safely neglected. If necessary, the condition 
Cm = 1 can be achieved by normalizing (y) over the 
observed frequency range (this measure would increase 
the noise and potentially introduce a bias, but it may be 
necessary to facilitate subsequent searches for a periodic 
signal) . 

The conversion probability -P 7 -> 7a is periodic in the 
wavelength domain, which requires that the model de- 
scription M(v) must not contain harmonic terms within 
a certain frequency range. This range is determined by 
several specific factors, including the distance to the ob- 
ject and the specific value of m 7a to be probed. This 
range is calculated and discussed below. 



A. Effective ranges of frequency and 
hidden-photon mass 

The sensitivity to detect the imprint of the amplitude 
a x with a signicance of n a on an observed radio spectrum 
with a Gaussian noise <r rms is given by a x = n a er rms (ne- 
glecting for the moment systematic errors resulting from 
an imperfect model representation, M.(y), of the astro- 
physical signal and systematic uncertainties of the mea- 
surement process). The highest radio frequency, v^, use- 
ful for recovering the oscillation signal can be estimated 
from Pr 



yielding (for n a > 1, using Eq. [2]) 



7T 

^ = 77 



2 arcsin(l/-i/ri~~) 



(12) 



1 /2 

with v\i = v± at n a = 1. At n a > 2, v^ pa (Tt/2)nJ v+ 
gives an estimate of v^ to within a 10% accuracy. 

The lowest frequency, v^ 1 containing a usable response 
from the oscillation is determined by the spectral spacing, 
Av. In this conservative estimate of vg is provided 

by the double of the Nyquist sampling rate, / s , of the os- 
cillation signal (this is required to take into account that 
measurements are made at a fixed set of frequencies and 
hence no "phase tuning" is feasible). The requirement 
corresponds to </? 7a (L, v t — Av/2) — <£> 7a (L, V£ + Av/2) = 
7t/2 and yields Vg = \J v±Av + (Av/2) 2 sa y/v+Av. 

Generic properties of a possible measurement are illus- 
trated in Fig. [TJ showing a modulation of an astrophys- 
ical signal and the frequencies vg and v^, as well as the 
observed frequency range 1^2] ■ 

The characteristic frequencies translate into an acces- 
sible mass range extending from m 7s to w 7a - The lower 
limit on the detectable hidden-photon mass m 7s is deter- 
mined by L = L osc (vi) (see Section III A[) . corresponding 
to v\ set to the frequency of the first local maximum 
^max.i- Tms yields 



m 



v Y \ 1/2 



MHz 



pc 



-1/2 



10-16 f^V 72 (L\ cV . (13) 



The largest accessible hidden-photon mass is ob- 
tained by setting vg = v%, which gives 



2tTV 2 



(LAv) 1 / 2 



1 v 2 f Av 
y/2v\ \v\ 



-1/2 



(14) 
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v,=(v,&.v) m 



1/A 




Hidden photon signal, 
F'(v) = F'(v)(l- 



log(v) 



FIG. 1. Modulation of an ideally modeled astrophysical sig- 
nal, F'(y) = 1, by the hidden photon oscillations. Measure- 
ments cover the [fiji^a] frequency range and are made with 
a spectral resolution of Av. Oscillations occur at a constant 
wavelength A*. Detection of the hidden-photon signal can be 
made if v\ < fmin,2 and 1/2 > za. Effective measurements 
(using the full sensitivity of the data) can be performed if 
v\ < vi and V2 > i>h- The spectral resolution must be better 
than v+jA. 



These limits are analyzed and presented in Fig. [2] for 
several existing and planned radioastronomical facilities, 
and for different types of astrophysical targets. In real 
experiments, the accessible mass ranges may be further 
limited by coherence effects and medium propagation 
as discussed in Sections III Al and IIIII This is illustrated 
in Fig. [5] by comparing the accessible ranges of hidden- 
photon masses to the limits imposed by the homogeneous 
ISM and IGM suppression. The impact of the medium 
suppression can be alleviated at lower photon masses by 
free propagation through a homogeneous medium (also 
illustrated in Fig. [2]), in addition to the favorable condi- 
tions that may exist for propagation through an inhomo- 
geneous medium. 

The sensitivity to x varies within the mass ranges de- 
scribed by Eqs. H3]fT4l The lowest hidden-photon mass 
m ^fuii £ or w j 1 j c j 1 a f measurements made in the 
[v\;v2\ frequency range is fully sensitive to x i s se t by 
the condition vg, = v±, which corresponds to 



J, full 



2-kvi 1 / Av 

(LAvf/ 2 = 71 W 



-1/2 



(15) 



The largest mass that can be detected with at the full 
sensitivity is then given by the condition v\ x = , result- 
ing in 



m^ fu11 = 2^27 



arcsm 



1 



1/2 



(16) 



with m^ fu11 = 2tt(ia j /£) 1/2 for n„ 
2\/27r(^2/ 'L) 1 ' 12 '"rhy 1 ^ for n„ > 2. With this approxima- 



1 and m"> fu11 

7s 



tion, m"^ ful1 can be expressed through m l 



, full 



V2 

,V4 I V X 



1/2 



(17) 



The mass limits described by Eqs. IT51[To1 can be used to 
estimate the frequency dependence (or, conversely, the 
hidden-photon mass dependence) of the sensitivity pro- 
vided by a given set of measurements for detecting the 
hidden-photon signal. 

Eqs. [15] and 1161 in particular, can be used for exper- 
iment optimization by requiring that m"^ > 
which corresponds to the inequality 



,£,ful] 



1/4 



1/2 



Av 



1/2 



> 1 



and connects the main parameters of the observational 
setup. For the goal of extending hidden photon studies 
to progressively lower hidden-photon masses, the most 
efficient strategy would therefore be to reduce v\. Im- 
proving the frequency spacing can result in producing a 
progressively larger number of datapoints that could not 
be used for probing higher hidden-photon masses. Hence 
the overall range of full sensitivity of a given experimen- 
tal setup would be reduced in this case. 



B. Sensitivity for the mixing angle \ 

The amplitude term a x of P^^ 1b implies that an n a 
bound on x ~ \/'V^rnis/2 can be obtained from multi- 
frequency flux density measurements. For N individ- 
ual flux density measurements with fractional errors a- u 



- = (££iof)~'7tf « °/N 1 ' 2 (if <n « a for all 
i = l,..., N). One can assume, as an example, that the 
astrophysical signal from a target object can be described 
by a simple power-law spectrum F(y) = F r (u/u T ) a and 
that the frequency dependence of the errors on flux den- 
sity measurements changes can also be described by a 
power-law dependence a{v) = a x {y jv^ (here, F r and a r 
refer to a flux density and its associated error, measured 
at an arbitrary reference frequency v T chosen inside the 
relevant frequency range). 

With these assumptions, measurements over the entire 
relevant frequency range [i/i>]i/ h ,} can be described by a 
characteristic signal-to-noise ratio 



1/2 



S = 



F{y) 



Av 



F T 



-it 1 -4 



a T vl (e+l)K'-^0 



(18) 

where £ = a — (3 and the integration limits are given by 
V(i = max(i/i,i^) and vy = min(^2,^h)- At a frequency 
spacing Av, the number of measurements contributing 
to the detection is iV meB = (iv — v^/Av, hence the 
effective cumulative signal-to-noise ratio of the data set 
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FIG. 2. Ranges of hidden-photon mass, m 7a , that can be probed with various radio astronomical instruments at different 
observing frequencies. The mass ranges are calculated for measurements done with planets (blue shades) , supernovae remnants 
(yellow shades) and active galactic nuclei (green shades). For each individual color, darker shades mark the mass ranges 
accessible to measurements with the first and second phases of Square Kilometer Array (SKA). The calculations are made 
assuming typical instrumental setups and generic ranges of distances to planets (0.5-10 au), supernova remnants (1-10 kpc) 
and active galaxies (0.02-3 Gpc). Brown dot-dashed lines mark the lower limits on detectable mass imposed by homogeneous 
ISM and IGM suppression. The dashed blue line illustrates the lower limit on detectable mass arising from propagation throuh 
an inhomgeneous IGM containing large scale voids (for the assumed void diameter of 100 Mpc). 



1/2 r 

is S = Nmes S. Recalling that, after accounting for the 
astrophysical signal, S~l/ cr Tms gives an estimate of the 
lowest achievable bound on the hidden photon coupling 



\lc 



<r t v$ (£+ 1)K- 



U/2 



(Av 



,1/2 



1/2 



4F r 



(19) 



The bound xiow remains constant for hidden photon 
masses m 7a fu11 < m 7a < m"^ fu11 and decreases rapidly 
outside this mass range as a progressively larger fraction 
of the measured data points are rendered outside the use- 
ful ranges of frequencies. 

Figure [3^ presents limits Xiow calculated for several 
existing and planned radio astronomical facilities and for 
measurements made with Mars (with F T = 30 Jy and 
a = 0.7 at v x = 86 GHz; |6l|), the supernova remnant 
Cassiopeia A (F r = 3000 Jy, cr rms = 100 Jy, v x = 1 GHz, 
a = -0.8; [H), and a fiducial compact AGN (F r = 10 Jy, 
a = —0.1, v r = 5 GHz) at a distance of 1 Gpc. These 
limits are obtained for the vacuum oscillation regime, 
without taking into account the potential medium sup- 



presion at lower hidden photon masses. For each of the 
instruments included in the plot, conservative assump- 
tions for generic technical parameters (summarized in 
Appendix B) have been made, hence the actual limits 
could be further improved by optimizing observations 
with a given telescope, for instance by increasing the 
spectral resolution of the measurements or applying ac- 
curate in-band (bandpass) calibration [63[. The figure 
clearly reflects the effect of improvements of cr rms (factor 
of ~ 10) and Av/v (factor of ~ 100) that will be pro- 
vided by the JVLA and SKA precursors, as well as the 
extension (factor of ~ 10) to lower frequencies provided 
by LOFAR. 



C. Source stacking 

Since the frequency behavior of the hidden-photon sig- 
nal is determined solely by the hidden-photon mass and 
the distance to the target object, signals from any num- 
ber of objects with known distances can be stacked to- 
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FIG. 3. Expected limits on \ that can be obtained with different radio astronomical facilities and different astrophysical objects 
under the assumption that the measurements are not affected by the plasma propagation effects (vacuum oscillations). The 
limits are calculated assuming generic instrumental parameters and: (a) single target sources at distances of 1 Gpc, 3.4 kpc 
[6(3 ] and 2 au for the AGN, Cas A, and Mars, respectively, (b) stacked data limits obtained from simulated populations of 100 
AGN (L =0.02-3 Gpc, F =1-30 Jy), 10 SNR (L =1-10 kpc, F =10-100 Jy) and 100 measurements made for the outer planets 
(L =0.5-10 au, F =10-200 Jy), with 7 = —2 used for all three populations. 



gcthcr to improve the resulting bound on the coupling 
constant. For the stacking, a suitable reference dis- 
tance L' can be chosen, yielding for each object located 
at a distance L a modification of observed frequency 
v' =v{L'/L). 

If all of the stacked sources have the same spectra and 
distances, stacking of iV bj objects will lead to a N^f 
improvement of xiow Let the stacked objects be drawn 
from a population with a uniform spatial density over dis- 
tances [I/ m i n ; -^maxL similar spectral indices, and an ob- 
served source count, N(F) = n(F r )(F/F r ) 7 over a range 
of flux densities [F min ; F max ]. 

In this case, stacking of iV Q bj spectra (each obtained 
with the same observational apparatus as described 
above) would yield a bound 

Xstack = XlowA^r 174 , (20) 

where F is estimated at a frequency v at which the signal- 
to-noise ratio S is achieved. Calculation of S may now 
involve different integration limits, as the frequencies vg 
and Vh must be calculated for a characteristic distance 
L. This distance is given by [(L 2 lax + L^ nin )/2] 1 / 2 for 
planets and Galactic objects and by [(if nax + L^ lin )/2] 1 ' 3 
for extragalactic objects. 

For the source and observation properties spec- 
ified by F{y) and a{v) (see Sect. IhTB|) , v = 
v^Sa./Frf/^, giving F = F^ a/i (Sa r ) a ^ . Con- 
sequently, N(F) = n(F T )(F / F T )i , where n(F T ) = 

N obi F?(l + 7) (F^ - F£?) ~\ 

Potential improvements of source stacking are illus- 
trated in Fig. which shows the limits on \ that can 
be achieved by stacking together measurements made for 
100 AGN, 10 SNR, and 100 measurements obtained for 



different outer planets of the Solar System. Improve- 
ments in both x an d the range of accessible hidden- 
photon mass are visible, compared to the single object 
limits in Fig The predictions for best cumulative 
limits from radio measurements in the entire 30 MHz to 
1400 GHz range are compared in Fig. 2] to the bounds 
derived from other experiments and observations. 

The potential effect of propagation through refractive 
media on the hidden-photon limits is illustrated in Fig- 
ure [5] for three different cases describing the minimum 
values of the plasma density, n . m ; n (Z), as a function of 
distance, I, along the propagation path. The three sce- 
narios for ?i c ,min(0 are adopted here solely to assess the 
potential range of possible outcomes of the propagation 
of the hidden-photon signal through refractive media. 

In all three scenarios, n eim i n (lkpc) = 10~ 5 cm~ 3 is 
adopted, implying effectively that the LOS paths for 
galactic objects cross at least one of the galactic "mini 
void" regions. The worst case scenario assumes that 
the photon signal from extragalactic objects propagates 
mostly through galactic ISM, i.e., not crossing large 
patches of IGM. In this case, n 0:m j n decreases from 
10~ 5 cm~ 3 at I = lkpc to 10~ 7 cm~ 3 at I = 3 Gpc im- 
plying that, for more distant objects, the LOS path has 
a progressively larger probability to cross lower density 
IGM regions. The average case corresponds to propaga- 
tion through typical IGM resulting in a n c . m ; n ~ 10 -5 - 
10~ 8 cm~ 3 range of densities. The best case allows for 
propagation through cosmic voids, with the correspond- 
ing n e ,min ~ 10~ 5 -10~ 10 cm~ 3 . For the planets, the den- 
sities of 50, 30, and 10 cm -3 (measured at 1 au and scal- 
ing oc (Z/lau) -2 ) have been adopted for the respective 
scenarios, based on measurements from [frij . 

The effects of resonant enhancement and medium sup- 
pression are clearly visible in the modified limits shown in 
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FIG. 4. Compound limits on \ expected to be achievable from observations at 0.03-1400 GHz compared with the limits obtained 
presently with other facilities and experiments 0, and references therein]. The limits from single object (blue) and multiple 
object stacking (red) are shown. Radio observations (particularly at frequencies below 40 GHz) will provide a unique probe for 
the hidden photon with masses below 10 _14 eV and extending down to ~ 10~ 17 eV where the measurements are likely to be 
limited by the medium suppression of the hidden-photon signal in the IGM plasma as described in sections III BHH Cl 



Fig. [5] Assuming that the best case scenarios would ap- 
ply for the majority of the lines of sight (since it is likely 
that a photon beam from a distant galaxy crosses one 
or more rarefied IGM or void regions), it is reasonable 
to conclude that the medium suppression would reduce 
the detectable \ to > 0.01 only for estimates made for 
m 7a < 10 _17 eV. Above this mass, the propagation ef- 
fects should not pose severe problems for constraining \ 
(and they indeed may even play a constructive role at 
least for some fraction of the photon mass range). 

D. Detection of periodic modulations 

The modulations induced by the hidden-photon signal 
on the broad-band spectrum of an astrophysical object 
can be best detected in the wavelength domain, where 
the modulation is sinusoidal, with the period given by 
A*. A fast Fourier transform (FFT) can be applied for 
searches in the data with dense and uniform coverage of 
wavelength space. For sparsely sampled data, epoch fold- 
ing 16511 or generic uniformity tests such as the Rayleigh 
test [66| or test [13] can be applied. 

These searches employ Eq. [11] with the Cm = 1 nor- 
malization of the residual flux density, which yields a 
functional form 

fx = /*[! + a* sin(u;*A + (f>o)] (21) 



of the periodic signal to be searched for. Assuming 
that the residual errors after the normalization are S\ = 
Crms /Cmi the parameters of the functional form are re- 
lated to the properties of the hidden-photon signal as 
follows: /* = 1 - a x /2 + 5\, a* = a x /(2/±), = 7Tf*, 
0o = 7r/2. Appendix [Al provides specific details of ap- 
plication of three different methods (FFT, epoch folding 
and generic uniformity tests) to searches for periodic sig- 
nals due to hidden photon oscillations. 

IV. DISCUSSION 

The analysis described above demonstrates principal 
feasibility of searching for a hidden-photon signal in 
broad-band radio spectra of cosmic radio sources, with 
galactic supernova remnants and radio-loud AGN pre- 
senting the best opportunity for extending the measure- 
ments below the hidden-photon mass of 10~ 14 eV, where 
essentially no measurements have been previously made. 

The mass ranges and kinetic mixing limits accessible 
for these potential searches are determined by several 
factors, including the oscillation and coherence lengths 
of the hidden-photon signal, the instrumental bandwidth 
and resolution, as well as the plasma density changes 
along the line-of-sight path to the source of the photons. 
The combined effect of the signal coherence and oscilla- 
tion length is expected to limit the searches to hidden- 



10 



-20 -18 -16 -14 -12 -10 -8 




-20 -18 -16 -14 -12 -10 



log 10 m yj [eV] 



FIG. 5. Modifications of the combined limits on \ due to propagation through refractive media. The dashed red line shows the 
limits predicted for object stacking observations under the vacuum condition (the same as the dashed red line in Fig. [4j. The 
best (solid black) , average (dashed black) and worst (dotted black) propagation scenarios illustrate potential effects of the ISM 
and IGM plasma for constraining the hidden-photon signal. For the best case scenario, the effect of resonant enhancement is 
visible. The medium suppression affects strongly the limits on % f° r hidden-photon masses below ~ 10" 17 eV. 



photon masses > 10 _19 eV, while the limits imposed by 
the propagation through the cosmic plasma may increase 
this limit up to ~ 10 _17 cV. 

Bandwidth and spectral resolution of the existing and 
planned radio astronomical facilities can support searches 
for hidden photons with masses well below 10~ 20 cV. The 
limiting instrumental aspect is the accuracy of ampli- 
tude calibration of radio receiver, which may limit plau- 
sible constraints on kinetic mixing to x ~ 10~ 3 . This 
problem may be alleviated by the advent of ultra broad- 
band receivers supporting in-band measurements across 
bandwidths in excess of 1 GHz. For this type of mea- 
surements, accurate bandpass calibration could deliver 
in-band amplitude accuracy of > 0.01%, thus potentially 
further lovering the limits on the kinetic mixing by at 
least an order of magnitude. 

For hidden-photon searches based on such in-band 
measurements, a combination of targets located at sig- 
nificantly different distances can be employed, profiting 
from the distance scaling of the hidden-photon signal. 
This ensures that any oscillatory pattern associated with 
a specific photon mass will be detectable only in one of 
the two measurements. This is realized for two objects 
at a distance ratio of L 2 /Li(L 2 > L\) > (n 2 /4)(n a /S v ), 
where 5 V is the fractional spectral resolution of the mea- 
surements. With such arrangements, the upper limit for 
an amplitude of a periodic oscillation in the bandpass 



obtained by dividing one of the two measured signals by 
the other would enable constraining x f° r the ranges of 
photon mass probed with either of the two targets. 

The effect of propagation on the sensitivity introduces 
a dependence on the assumed electron density n c of the 
intervening medium and therefore on the line of sight. 
For galactic sources, estimates of n e at a specific line of 
sight can be obtained from pulsar dispersion measure- 
ments [52l , [53! ] . Optical hydrogen absorption lines can be 
used for assessing the line-of-sight structure of the IGM 
[68j , which can be applied for measurements in individual 
extragalactic targets. The situation may improve sub- 
stantially after the large-scale HI surveys planned at the 
LOFAR, MeerKAT, ASKAP and the SKA (gUzl would 
deliver a very detailed picture of the IGM up to very high 
redshifts. 

In the absence of information about the plasma prop- 
erties on individual lines of sight, a potential remedy for 
the object stacking would be to make generic assump- 
tions on the minimum plasma density as a function of 
distances to individual targets (similarly to the approach 
employed in Section [ill Cp . Despite being inherently im- 
precise, this approach should still enable improving the 
constraints obtained on x from object stacking. 

Data from the next generation of large scale continuum 
and spectral line surveys at radio wavelengths will pro- 
vide sufficiently accurate information about the broad- 
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band continuum, line-of-sight distribution of the IGM, 
and distances to many thousands of extragalactic radio 
sources. This will make object stacking the tool of choice 
for the hidden-photon searches in the radio regime and 
will certainly lead to strong improvements of the limits 
obtained from radio data on the kinetic mixing of hidden 
photons in the 10~ 14 -10~ 17 cV mass range. 
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where tfi = 7ri/(2n D ) and £ w = 0.773 [Hf is a correction 
factor taking into account the finite bin width. Reaching, 
in a given bin, a desired confidence level, c (in percent), 
of detection of the oscillations, implies (wi) > Wq, where 
Wo is derived from the % 2 probability distribution, P2(x 2 ), 
with two degrees of freedom, requiring that 

/>oo 

1 - c/100 = (A/AO / P2(xW- (A5) 

If a signal is detected, the amplitude a* (and conse- 
quently, the kinetic mixing parameter a x ) can be ob- 
tained by requiring that 

poo 

c/100= / P2 (xW, (A6) 
with w a = wq — (wt) + 8/2. 



Appendix A: Analysis techniques for periodic signals 

1. FFT searches 

An FFT search can be employed effectively for recover- 
ing the hidden photon signal, if measurements are made 
in a frequency range [1/15^2], corresponding to a wave- 
length range A = Ai — A2. The measurements are as- 
sumed to be sampled densely enough to be binned into 
?ib = 2 m bins, with each bin described by the bin flux 
density /, and its associated error 5i. This situation can 
be realized, for instance, for observations with high spec- 
tral resolution (e.g., with LOFAR, JVLA, or ALMA) ap- 
plied to search for hidden photons with m 7s > 2iVy/ V2/L, 
for which > V2- The resulting 2m — 1 Fourier coeffi- 
cients are given by 



E 



fj exp(iui\j) 



(Al) 



where Xj = (j — l/2)A/rib is the arithmetic mean wave- 
length of the j bin and fj is the respective bin-averaged 
flux density. The corresponding power spectrum Wi is 
given by 



w. 



2\a t \ 2 /F u 



(A2) 



where Ft, 



"tot = X)j=i fj * s ^ e total flux density in the 
bins. The average power spectrum contribution from the 
measurement noise is accounted for by the term 



~s = -ys, 



(A3) 



Adding the noise contribution and taking into account 
frequency dependence of \a,i\ 2 [65|, the average power in 
the i th bin can be written as 



a* sin (fii 



(A4) 



2. Epoch folding 

If the data arc not homogeneously sampled across 
the measured wavelength domain A but still can be 
divided into phase bins of size < Av jvi (thus giving 
TjU > 2^2(^2 — v\)l(y\ Av) phase bins), epoch folding [c/., 
lola . F73l ] can be effectively performed in order to search for 
all hidden-photon masses in the [m*;m"J range. After 
epoch folding with a given trial wavelength A', each bin 
is characterized by the bin flux density fa and its asso- 
ciated error 8i . The significance of the signal can be as- 
sessed using the statistics s = Ym=i (/« — f) 2 /^h where 
/ = Ftot/ritp. For the same assumptions as used above 
for the FFT searches, the mean value of the statistics is 
given by 



a* sin ip 



(A7) 



with ip = 7r/?i0 and £0 = 0.935 [65|]. The sensitivity 
to oscillations is established at a confidence level c by 
requiring (s) > sq where sq is obtained from the condition 



1 - c/100 = (A/A') / Pn,-i(x 2 )d X 



(A8) 



Similarly, the amplitude of the detected signal can be 
estimated from 



c/100 = / p n ,_! (xW, 



(A9) 



where s a = sq — (s) +6/2. 



3. Generic uniformity tests 

For poorly and unevenly sampled data, generic unifor- 
mity tests such as the Rayleigh test [§6|, [73, [74[ , the 
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TABLE IV. Generic technical parameters of radio telescopes 
used in the calculations of sensitivity for \ m the radio regime. 



Telescope 




V2 


iVch 




p 




[GHz] 


[GHz] 




[Jy] 




LOFAR 


0.03 


0.24 


1000 


0.3 


-0.25 


SKAi 


0.08 


0.33 


5000 


0.04 


-0.1 


ASKAP 


0.1 


1.4 


2000 


0.1 


-0.25 


SKA 2 


0.3 


3.0 


5000 


0.04 


0.1 


MeerKAT 


0.6 


10 


200 


0.03 


0.25 


Effelsberg 


0.3 


37 


40 


0.1 


0.25 


JVLA 


0.3 


13 


500 


0.1 


0.25 


APEX 


170 


410 


4000 


10 


0.25 


ALMA 


86 


1389 


5000 


1 


0.25 



moments entering the statistics 



f'i cos(j (/>i) 



^//sin(i^) 



(All) 

The resulting calculated powers 2nR 2 and Z^ n should be 
tested against \ 2 distributions P2(x 2 ) an d P2m(x 2 )- rG_ 
spectively. Adopting the same approach as for the epoch 
folding, these values will yield confidence limits for de- 
tecting a periodic signal with the wavelength A'. 



Appendix B: Basic technical characteristics of 
simulated radio observations 



test [67J or the H-test [75| can be applied, relieving also 
the requirement to bin the data before searching for a 
periodic signal. 

Similar to the case of epoch folding, the data con- 
sisting of n flux density measurements have to be first 
normalized by a factor of / m i n = rnin(/j) and folded 
with a trial wavelength A', yielding a set of amplitudes 
f'i = fill 'min and phases fa, with i = 1, .., n. Calcula- 
tion of the respective trigonometric moments is done us- 
ing the terms f- cos fa and sin fa, so that the resulting 
Raylcigh power becomes 



n R z = - 

n 



COS ( 



(A10) 

Similar modification should be done to the trigonometric 



Tabic IIVI describes general technical parameters 
adopted for simulating observations with the radio tele- 
scopes used for the predictions of the hidden-photon mass 
ranges and the limits on \ from measurements in the 
radio regime at frequencies of 0.03-1400 GHz. The pa- 
rameters presented in the table are: the lowest v\ and 
highest V2 observing frequencies, the number of measure- 
ments, iV c h, available withing the observing range, the 
r.m.s. noise, cr r of a single measurement at the reference 
frequency v r , and the power index (3 describing the fre- 
quency dependence of the r.m.s. noise. The parameter 
values given in Table IIVI provide only basic benchmark 
description of technical capabilities of the individual tele- 
scopes. Conservative estimates have been adopted for the 
number of spectral channels and the r.m.s. noise (assum- 
ing typical "shallow" survey observations), and these val- 
ues can be improved by one or more orders of magnitude 
by employing full correlator capabilities and dedicated 
targeted observations. 
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